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Abstract 



Hadron-nucleus collisions at LHC energies are studied by including explicitly 



semi-hard parton rescatterings in the collision dynamics. Under rather general con- 
£T} ■ ditions, we obtain explicit formulae for the semi-hard cross-section and the inclusive 

minijet transverse spectrum. As an effect of the rescatterings the spectrum is low- 
ered at small pt and is enhanced at relatively large transverse momenta, the defor- 
mation being more pronounced at increasing rapidity. Its study allows to test the 
proposed interaction mechanisms and represents an important baseline to examine 
Cu nucleus-nucleus collisions. 
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1 Introduction 



Given the rapid growth of the hard cross-section in hadronic and nuclear collisions [||] , the 
typical inelastic event will be dominated by the perturbative regime at very high energies 
so that, at the LHC, one may expect to be able to derive global features of the inelastic 
interaction by perturbative methods. Such a capability, unavoidably limited to a restricted 
number of physical observables, implies however a few non trivial improvements in the 
understanding of the mechanisms operating in the interaction process. To be estimated 
in a sensible way, different physical quantities may in fact need a different degree of 
understanding of the interaction dynamics, since many details of the process may be of 
little relevance for some observables, while they may be essential for other quantities. 
Identifying and evaluating such physical observables represents a non trivial improvement 
in our capability of using the perturbative QCD to describe physical processes. 

An obvious problem will appear when trying to elaborate along these lines. A pertur- 
bative calculation does not introduce any scale in the dynamics, so that in this case the 
kinematic variables are the quantities which give the dimensionality to the related physical 
observables. The property is associated to the relatively low rate of events falling in the 
domain of the canonical fixed-x perturbative QCD. On the other hand the dimensional 
factor which characterizes the global features of the typical inelastic event is, rather, the 
hadron (or nuclear) scale. When the perturbative regime dominates a physical observable 
which represents global features of the inelastic interaction, the hadron (or nuclear) scale 
should therefore appear also in the corresponding perturbative calculation, presumably 
introduced through the non-perturbative input. The structure functions, namely the up 
to now non-perturbative input of basically all perturbative calculations, are on the other 
hand dimensionless quantities. This implies that the structure functions, in their present 
form, will no longer be an adequate non-perturbative input when trying to accomplish 
the program outlined above. 

A related aspect is the complexity of the interacting states. The canonical, fixed- a;, 
perturbative QCD approach considers only perturbative processes initiated by a pair of 
partons. The approach is appropriate in the case of very dilute interacting systems, while 
it becomes obviously inadequate in a regime with very large parton densities. In the 
case of a partonic interaction in the black disk limit, the initial configuration is in fact 
isotropic in transverse space, differently from the final state produced by an interaction 
initiated by two partons (namely, at the leading order in as, two jets back-to-back in 
Pt), where a direction in the transverse plane is singled out. A natural way to recover 
the black disk symmetry in the final state, is to include in the interaction perturbative 
processes initiated by more than two partons (namely, semi-hard parton rescatterings) , 
whose relevant property is to produce many large-p t jets also at the lowest order in a$- 
A non-trivial feature is the associated non-perturbative input. To deal with processes 
initiated by more than two partons one needs in fact to introduce, as a non-perturbative 
input, the many-body structure functions, which contain independent informations on the 
hadron (or nuclear) structure with respect to the one-body structure functions needed to 
deal with processes initiated by two partons. A basic property is that the n-body structure 
functions are dimensional quantities, in such a way that when n is larger than one the 
many-parton initiated processes introduce non-perturbative scale factors in the dynamics 
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of the interaction in a natural way, allowing one to deal with the problem of dimensionality 
previously mentioned. 

By introducing interactions initiated by may partons one may therefore gain the capability 
of describing, by means of perturbative QCD, at least a few general properties of the 
typical interaction at the energy of the LHC. To pursue such a program one should then 
i) evaluate in perturbative QCD processes involving many partons in the initial state, ii) 
face the problem of the unknown non perturbative input and develop a strategy in that 
respect, and Hi) study the infrared problem by finding observable quantities which are 
infrared stable. This last step represents the final achievement of the whole program. 

The purpose of the present paper is to discuss the case of hadron-nucleus interactions 
(hA, for brevity). Being intermediate between hadron-hadron (hh) and nucleus- nucleus 
(AA), hadron-nucleus interactions allow several simplifications in the formalism developed 
to discuss heavy- ion collisions. In fact, as it will be shown hereafter and differently with 
respect to the latter case, in the hadron-nucleus instance we will able to obtain closed 
analytic expressions for the semi-hard cross-section under rather general conditions. We 
will then study the inclusive minijet transverse spectrum, which is related in a direct 
way to the underlying dynamics and is therefore an important baseline for the study of 
nucleus-nucleus collisions. 

Beside its intrinsic interest, inclusion of semi-hard rescatterings in the computation of the 
transverse spectrum has been advocated by many authors [| f|, |^, |5| as the basic mech- 
anism underlying the Cronin effect 0, namely the deformation of the hadron pt-spectra 
in nuclear collisions as compared with the expectations of a single large-p 4 production 
mechanism. Multiple parton collisions have also been related to higher-twist parton dis- 
tributions |§ H |H]]. A non-perturbative study of the transverse spectrum in hA collisions 
in the framework of the McLerran-Venugopalan model for nuclear and hadronic collisions 



was presented in [11 



Another reason of interest in hadron-nucleus collisions is that theoretical models can be 
tested against experimental data in a situation where further nuclear effects are absent, 
like, e.g., the formation of a hot and dense medium which can further modify the trans- 



verse spectrum via energy-loss [12, 13]. Therefore a detailed understanding of hA collisions 



represents an important baseline for the generalization to AA collisions [[14], |T5[ and for 



the discovery of novel physical effects [16 



An explicit approach to semi-hard interactions in heavy ion collisions at the LHC on 
the lines previously described, has been accomplished, at least partially, with the help 
of a few simplifying hypotheses. The program has been implemented in [|TT| , |18| , [19| , |20| , 
and various physical quantities have been evaluated in |H], p2| . The approach relies, to 
a large extent, on the idea of self-shadowing, which we recall for completeness in Sec. |2|. 
In Sec. ^| we discuss our expression of the semi-hard interaction probability between two 
colliding partonic configurations. We discuss also the multi-parton distributions, that are 
studied with a functional formalism and finally combined with the interaction probability 
to derive the hadron-nucleus semi-hard cross-section. Sec. |] is devoted to the discussion 
of the inclusive minijet transverse spectrum, with particular emphasis on the mechanism 
of subtraction of infrared divergences, which is explicitly implemented in our approach. 
Results of numerical evaluations of the inclusive spectra of minijets in hadron-nucleus 
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collisions are presented in Sec. |5]. The last section is devoted to the concluding summary. 



2 Self-shadowing 



To face the problem of unitarity corrections we make use of the self-shadowing property 
of the hard component of the interaction. For the sake of completeness, in the present 
paragraph we recall the main points about the self-shadowing cross-sections in hadron- 
nucleus interactions 



Let's consider the inelastic hadron-nucleus cross-section {(Tin) a, whose expression may 
be expanded, in the Glauber approach, as a binomial probability distribution of inelastic 
nucleon-nucleon collisions: 



[O'inJA 



d 2 /?[l - (l- 
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(2.1) 



In Eq. (|2.1| ) r(j3) is the nuclear thickness function, which depends on the impact parameter 
(3 and is normalized to one, A is the atomic mass number and o in is the inelastic hadron- 
nucleon cross-section. One may classify all events according to a given selection criterion, 
which we call C, while we call TV the events that are not of kind C. We assume that in a 
hadron-nucleon collision all events of kind C contribute to ae, all other events contribute 
to <7/v, so that the inelastic hadron-nucleon cross-section may be written as 



Or 



One may then ask for the expression of the cross-section (ctc)a to produce events of kind 
C in a collision of a hadron against a nuclear target. Then, to obtain (<Jc)a, one may 
express (<Ji n ) n in Eq. fl2.1|) as a binomial sum of "elementary" events of kind C and of kind 
Af: 



k=0 ^ ' 



(2.2) 



An interesting case to consider is when the events of kind C are such that any superposition 
of elementary events of kind C, both with events of kind C and of kind J\f, always gives an 
event of kind C. In this case, all the terms of the sum in Eq. ( |2.2|) , with the only exception 
of the term with k = 0, contribute to (ctc)a, which is therefore given by: 



[vc)a 



n=l ^ ' Lfc=l V / 



By using the relation 
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one obtains: 



A 



A-n 



n=l ^ ' 

d 2 f3\[a in r(f3) + 1 - a m r(f3)Y - (a^r{P) + 1 - ^{(5)) 
d 2 p[l-(l- aerify 



n=l 



A 



(2.3) 



Notice that, in spite of the fact that we included superpositions of elementary events of 
kind C with events both of kind C and of kind Af, the nuclear cross-section (ctc)a is obtained 
by summing all possible multiple hadron-nucleon interactions of kind C alone with a 
binomial probability distribution, precisely as {(T in )A is obtained by a binomial distribution 
of hadron-nucleon inelastic interactions. This relation states the self shadowing property 
of the events of kind C: all unitarity corrections, namely the term [l — oct(/3)] in the 
third line of Eq. (|2.3| ), are expressed by means of the cross-section oc only. However, 
this does not mean that [pc)A doesn't contain events of kind Af, but rather that they 
are irrelevant to obtain (ctc)a- The property that an event of kind C remains of kind C 
even after any number of events of kind Af translates into the disappearance of in the 
nuclear cross-section ((?c)a- 

Given the discussion above, the only part of the nuclear interaction that still misses is 
the cross-section for elementary events of kind Af alone. It can be obtained by considering 
the following difference 



d(a. 



in) A 



d 2 (3 



d 2 (3 



1 - a c r{f3) 



1 - (t C t(P) 



i - (<T C + wMP) 



X <1 



1 - cr c r((3) 



(2.4) 




A-k 



1 - acr(P) 



which is therefore bounded by [l — o"c r (/5)] (second line of |2.4| ), namely by the proba- 
bility of not having any interaction of kind C at a given impact parameter j3. The ratio 
0jv"7"(/3)/[1 — (Jct((3)] is in fact a quantity smaller than one, since a in r(P), which is equal 
to (cr c + tj^)t(P), is a probability. It may be understood as the probability of an hadron- 
nucleon interaction at a given impact parameter, under the condition that no event of 
kind C takes place. Hence the last line of Eq. ( |2.4j) shows that after removing all events 
of kind C the interaction is expressed by a binomial distribution of events of kind Af. 

Finally, we observe that if we compute the average number of hadron-nucleon collisions 
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of kind C, (n)(ae)A, rather than the cross-section (<tc)a, the result is: 
(n)M A = j d 2 (3^n( A \{a C T{P))\l-a c r{(3)) A - n 

J n=l ^ ' 



7=1 



Notice that the average number of interactions of kind C is expressed by the single- 
scattering term, without any unitarity correction. 



3 Semi-hard cross-section 



In this section we want to represent the semi-hard hadron-nucleus cross-section analo- 
gously to the self-shadowing cross-section (|2.3| ), but considering as elementary objects the 
partons instead of the nucleons. Indeed, the hard component of the interaction satisfies 
the requirements of the self-shadowing cross-sections if one assumes that a parton which 
has undergone interactions with large momentum exchange can always be recognized in 
the final state. An immediate difference with respect to the previous case is that now 
there is no upper bound on the number of partons that can take part in the collision. 
Because of self-shadowing all unitarity corrections to the semi-hard cross-section will be 
therefore expressed by means of the semi-hard partonic cross-section only, so that one 
doesn't need to make any commitment on the soft component when only the semi-hard 
part of the interaction is of interest. Self-shadowing allows moreover to control also the 
soft component of the interaction by perturbative means, since that contribution is lim- 
ited to a fraction of the cross-section proportional to the probability of not having any 
hard interaction at all (see Eq. |2.4| ). Obviously the unavoidable restriction of all consider- 
ations done by perturbative means is that those are limited to partonic final states, whose 
properties will hopefully survive hadronization. 

To represent the interaction between hadrons and nuclei in terms of partonic interac- 
tions, each one with relatively large momentum exchange, one needs to write the cross- 
section for a given non-perturbative input, namely for a definite partonic configuration 
of the two interacting objects. Then, as a perturbative input, one needs to write the 
probability of having at least one semi-hard interactions between the two configurations 
of partons. We discuss these two inputs in the next two subsections and in Sec. |0] we 
combine them to obtain the hadron-nucleus cross-section. 



3.1 Perturbative input: semi-hard rescatterings 

In a processes involving many partons in the initial state, an important distinction is 
between connected and disconnected hard interactions. A hard process may in fact be 
represented also by a disconnected hard amplitude, since the overall interaction process 
may still be connected by the soft component of the interaction. The simplest case of a 
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disconnected process is obviously the one where the hard part is represented by two par- 
tonic interactions at the lowest order in perturbative QCD, corresponding to two 2^2 
parton scatterings. Since all hard collisions are characterized by short transverse dis- 
tances, disconnected hard processes give rise to a picture of the interaction where the 
different partonic collisions are all localized in different points in the transverse space, 
with a transverse distance of the order of the scale of soft interactions ||24j| . The discon- 
nected component of the hard-interaction leads therefore to a geometrical picture of the 
process, giving as well some indications on the degrees of freedom characterizing the non 
perturbative input. The many-body parton distributions need in fact to depend explic- 
itly on the parton transverse coordinates, to identify the partons involved in each given 
sub-processes with a definite localization in transverse space. Notice that the dependence 
on the transverse coordinates and the number of partons taking part the interaction are 
the basic information needed to assign the dimensionality to the many-body structure 
functions, and therefore to introduce the non-perturbative scale factors in the interaction 
dynamics. 

While the main feature of the disconnected component of the hard amplitude is to 
give rise to a geometrical picture of the semi-hard interaction, the connected component 
of the amplitude becomes more and more structured when one approaches the black disk 
limit, where a single projectile parton may interact with several target partons with large 
momentum exchange in different directions in transverse space. The simplest possibility of 



such an interaction was discussed in Ref. [20|, where the forward amplitude of the process 



and all the cuts were derived in the case of a point-like projectile against two point-like 
targets, in the limit of an infinite number of colors and for t/s —>■ 0. In this case one finds 
that the different cuts of the 3 — > 3 forward amplitude are all proportional one to another 
and the proportionality factors are the AGK weights [p^[] . A consequence is that one may 
express the three-body interaction as a product of two-body interaction probabilities. The 
results obtained in that simple case may indicate a convenient approximation of the many- 
parton interaction probability. One can in fact argue that the many-parton interaction 
process may be approximated by a product of two-parton interactions, so that one can 
call the process re-interaction or rescattering. The whole interaction is therefore expressed 
in terms of two-body interaction probabilities, precisely as the interaction between two 
nuclei is expressed in terms of nucleon-nucleon collisions. Hence, given a configuration 
with n partons of the projectile and I partons of the target, we introduce the probability, 
V n i, of having at least one partonic collision, in a way analogous to the expression of the 



inelastic nucleus-nucleus cross-section ||26 

V n , = 



n IB 1 -*«)]' w 

i=l j=l 



where is the probability of interaction of a given pair of partons i and j. Since the 
distance over which the hard interactions are localized is much smaller than the soft in- 
teraction scale, one may approximate a(xiXf bi — bj) ~ a(xiXj)S ( - 2 \b i — bj), where Xi and 
Xj are the momentum fractions of the colliding partons, bi and bj their transverse coor- 
dinates and o~(xiXj) is the partonic cross section, whose infrared divergence is cured by 
introducing a regulator p . For example, po may be the lower cutoff on the momentum 
exchange in each partonic collision, or a small mass introduced in the transverse prop- 



7 



agator to prevent the divergence of the cross-section at zero momentum exchange. The 
expression of V n ,\ is the analogue of Eq. ( |2.3|) and represents the explicit implementation 
of self-shadowing for the interaction of two partonic configurations. 



3.2 Non-perturbative input: multi-parton distributions 

In this section we discuss the non perturbative input to the process. To approach the 
problem in the most general form we use the functional formalism introduced in ||19|| . At 
a given resolution, provided by the regulator p , one may find the nuclear (or hadronic) 
system in various partonic configurations. We call P^ 1 (u\ . . . u n ) the probability of a 
configuration with n-partons (the exclusive n-parton distribution) where Ui = (bi, Xi) 
represents the transverse coordinate of th i-th parton, fej, and its longitudinal fractional 
momentum, Xj. The distributions are symmetric in the variables u^, and can be obtained 
from a generating functional defined with the help of auxiliary functions J{u) as follows: 

Z\J] = S ^-] I J{u 1 )...J{u n )P {n) (u 1 ,...,u n )du 1 ...du n , 

all infrared divergences are regularized by po, which is implicit in all equations. Probability 
conservation yields the normalization condition Z[l] = 1. Then, the exclusive n-parton 
distributions can be obtained by differentiating the generating functional Z with respect 
to the auxiliary functions: 



P (n) («!,...,«„) 



8J{ Ul ) "' 5J{u n ) L J K=° 



A useful representation of Z may be found by introducing its logarithm, J 7 , with normal- 
ization JF[1] = 0, so that 

Z[J] = , 

and by studying the inclusive n-parton distribution, D^ n \ They can be obtained as func- 
tional derivatives of Z or of T . Indeed 



u 



P x {u) 

5Z 
6J(u) 



P( 2 \u,u')du' 
8F 



P (3 \u,u',u")du'du" + 
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Ui,U 2 ) 



P (2) (mi,m 2 ) + J P^ 3 \u 1 ,u 2 ,u' )du'+ l - J P^\ Ul ,u 2 ,u', u")du'du" . 



5 2 Z 



5J(u\)5J(u 2 ) 



b 2 T 



j=i 



SJ(ui)SJ(u 2 ) 



j=i 



SJ(ui) 5J(u 2 ) 



j=i 



and so on for higher multi-parton distributions. These relations show that the correlated 
part, C^ n \ of the inclusive n-parton distribution (also called n-parton correlation) is simply 
given by differentiation of the generating functional T: 



C (n \ Ul ,...,u n ) 



S 



8J{u\ 



8J{u n ) m \^ 
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so that the expansion of T near J = 1 reads: 



;F[J] = / I»[J(u) - l\du 

OO „ 

n=2 ' "* 

where T(u) = D^(u) for consistency with the notation used in previous papers. In this 
way we have obtained a convenient representation of the generating functional Z = exp[jF] 
in terms of the single parton inclusive distribution, T, and of the multi-parton correlations, 
C( n \ In the simplest case where we neglect all the correlations between the partons, 
namely C^™- 2 -* = 0, the generating functional is given by 

Z[J\ = e /r(«)[J(«)-i]*i _ (3.2) 



3.3 Hadron-nucleus cross-section 

The general expression of the semi-hard cross-section at fixed impact parameter is ob- 
tained by folding the interaction probability, Eq. (|3.1|) , with the multi-parton exclusive 
distributions of the two colliding systems (in our case a hadron, h, and a nucleus of atomic 
number A): 



d 2 P ~ J ^ n! 5J(m -(3)"' SJ(u n -(3) h[ 1 



j=o 



oo 

X 

m 



1 5 6 z fr/ll 

^ ml 5 J'K) " ' SJ'«) A[ J ^'=o 

n m n m 

x ji-nni 1 - & «oi } n ^ n ^ > ( 3 - 3 ) 

i=l J=l i=l j=l 

where (3 is the impact parameter between h and A. To simplify the notation we introduce 
the following operators: 

5j = [ dui—— — ; 8'a — [ dui 



5J( Ui -(3) ' 3 J dJ'iu'j) ' 

Given two functions / = f{u) and g = g(u), the following identity holds: 

e s -fZ[J + g] = Z[J + g + f] , (3.4) 

where 5 ■ f = f du5 /8J{u)f{u). In other words, the exponential of the operator 5 acts on 
the generating functional Z by shifting its argument of the amount /. 

In the case of hadron-nucleus interactions one may be allowed to neglect the rescat- 
terings of the partons of the nucleus. Indeed, even at very high center of mass energies 
the average number of scattering per incoming parton is smaller than the average number 
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of nucleons along the parton trajectory, except in the very forward rapidity region J2T . 
With this assumption the interaction probability may be simplified as follows: 



n,m 



{l - II [1 - Vij] } - 53 & V J\ Y Y &i ^ kl + ^ Y Y & H & klVrs + ■■■ (3.5) 
i,j i,j i,k j^l i,k,r jj^lj^s 

After contracting Eq. ( |3.5|) with the differentiation operators in Eq. ( |3.3|) one obtains 

q 



q>l 



i=l 



E—8i ...5 n 
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1 — exp^— > <) • rr, 



1 — exp 



r 



J+S' 



By using the identity 

d 2 p ~ 



the semi-hard cross-section becomes: 



1 — exp 



5 ■ (e~ 5 '- a - 1 



Z h [J+l}Z A [J' + l] 



(3.6) 



J=J'=0 



This result is very general and includes all possible parton correlations of both the pro- 
jectile and the target; the only assumption made is that target partons do not suffer 



any semi-hard rescattering (we will comment more on this assumption in Sec. |5.1|) . A 
meaningful approximation (see Ref. ||19||) is to consider the nuclear partons uncorrelated, 
namely C^~ 2 ' ) = 0. Then, by using Eq. ( p.2| ) the cross-section reduces to: 



da 



d 2 (3 



f =ESE(-ir m (l)^[i--]^+i]|^ 

n=l ' m=0 V 7 

1 - z h \e~f H '' u ' )rA{u ' )du '~ 



(3.7) 



If we neglect also the correlations between the partons of the projectile, we get a further 
simplification: 



dan 
~d?f3 



1 — exp< — / duTh(u — (3) 



I _ g-J a(u,u')r A (u')du' 



(3.8) 



Both in Eq. ( |3.7|) and in Eq. (|3.8j ) the cross-section is a function of 



W h {u,P) =T h {u-f3) 



1-e 



- / &(u,u')T A (u')du' 



T h (u-P)V A (u) 



(3.9) 



which represents the number of projectile partons that have interacted with the target, 
i.e., the projectile wounded partons |17|, we call them minijets, even if they did not yet 
hadronize. Va{u) represents the probability that a projectile parton with given u = (x, b) 
has at least one semi-hard interaction with the target, hence the cross-section is obtained 
by summing all events with at least one interaction. 



One might obtain the average number of wounded partons, Eq. ( |3.9| ) by working out 
directly from Eq. (|3.3| ) the average number of projectile partons which have undergone 
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hard interactions |17 , [19|| (a detailed numerical study of this quantity in Pb-Pb colli- 
sions at LHC and RHIC energies is presented in |21 , |22fl). The result, Eq. ( p.9|) , is 



obtained under the only assumption that all the target partons are uncorrelated. There- 
fore, f duW(u,/3) = (n)da/d 2 p represents the integrated inclusive cross-section to detect 
all scattered projectile partons, and takes into account the correlations of the projectile 
partons at all orders. Of course, the projectile parton correlations appear explicitly in 
the total hadron-nucleus cross- sect ion. In the simplest case of two-parton correlations one 
would obtain: 

^| = l-exp|-y duW h {u-i3) + ^ J dudu'V A {u)C^\u- f3,u' - 0)V A (u')} . 

(3.10) 

The effect of correlations on donjd 2 fi is however small, both when unitarity corrections 
are small (i.e., when the semi-hard parton-parton cross-section is small, so that Va an d 
Wh are both of order cr#) and when they are large (i.e., when an is large, Va ~ 1 and 
Wh is large). If, on the other hand, one is looking for correlations, the simplest quantity 

(2) 

which depends linearly on is the double-jet inclusive cross-section. 



4 Inclusive minijet transverse spectrum 

The whole semi-hard hadron-nucleus cross-section results from the superposition of the 
multiple interactions of the partons of the projectile hadron (which is a dilute partonic 
system) with the nuclear target (which is a dense partonic system). The inclusive trans- 
verse spectrum of the projectile minijets is given by the distribution in p t of the average 
number of wounded partons, and is affected by the presence of semi-hard rescatterings 
rj|. The deformation of the high-p t hadron spectra, which leads to the Cronin effect, 
was studied in terms of semi-hard parton rescatterings in |fj, ||, ||, [J, where partons that 
suffered up to two scatterings where included, leading to a good description of the data for 
pA collisions up to \fs = 39 GeV/A. However, the two-scattering approximation breaks 
down at higher energies, except at very high p t , and the whole wounded parton transverse 
spectrum is needed. More phenomenological approaches || |T1| , which take into account 
also intrinsic transverse momentum, model the effects of multiple scattering as an addi- 
tional Gaussian p r broadening for each rescattering suffered by a parton. A random-walk 



model of the multiple scatterings was proposed in |L5 . 

After the introduction of semi-hard parton rescatterings, integrated quantities like 
the semi-hard cross-section and the minijet multiplicity show a weak dependence on the 
infrared cutoff needed to regularize the infrared divergences arising in the perturbative 



computations [|17|, pi| . On the contrary, it will be shown that differential quantities like 
the minijet pt-spectrum are more sensitive on the detailed dynamics of the interaction and 
show a stronger dependence on the cutoff, if only logarithmic. To reduce this dependence 
on the cutoff one needs to improve further the picture of the dynamics by including also 
gluon radiation in the interaction process. Some steps along this line in the case of deep 



inelastic electron- nucleus scattering have been presented in ||10|| . In this paper, however, 
we neglect the problem of the gluon radiation and we concentrate on the effects of elastic 
rescatterings. 
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4.1 Transverse spectrum 



We can expand the average number of projectile wounded partons, Eq. (|3.9| ), at a given 
x and b in a collision with impact parameter (3, in the following way: 

W h (x, b, (3) = T h (x, b-(3)J2 {Mx ; b)r e-<-^> , (4.1) 



u=l 



where (ua(x, b)) = J cIx'TaIx', b)a(xx') is the average number of scatterings of a projectile 
parton at a given x and b ||19|| . The average number of wounded partons is then given 
by the average number of incoming partons, Th, multiplied by the probability of having 
at least one semi-hard scattering, which is given by a Poisson distribution in the number 
of scatterings, z/, with average number (ha{x, b)). Therefore, we can obtain the inclusive 
differential distribution in p t by introducing a constraint in the transverse momentum 
integrals that give the integrated parton-parton cross sections in the expression above: 

^(x,b,fl =r fc fo&-/3)f;± jv A {x^b)..S A {x' y Me- idx ' VA{x '^^ 
da da 



x ^-■■■^-S^(k 1 + --- + K-p t )d 2 k 1 ...d 2 k v dx' 1 ...dx' v . (4.2) 

d l k\ d A k v 

The limits of integration on x\ and x' are respectively xx[s > 4k 2 and xx's > 4pg, and all 
the distribution functions are evaluated for simplicity at a fixed scale. 

By using the above formula one can study the pt-broadening of a wounded parton, 
namely, the square root of the average transverse momentum squared acquired through 
its path across the nucleus. Consider a single projectile parton with fixed x and b. The 
probability that it acquires a certain p t after the collision is given by Eq. (|4.2|) divided by 
the number, Th(x, b — (3), of incoming partons: 



d?p t d?p t ' T h (x,b-/3) 
Then, the average transverse momentum squared of a wounded parton is given by 



(p t 2 (x,&)), 



«P?» 



((1)) ' 

where ((f(pt))) = J d 2 p t f{pt)dV 'a/ 1 d 2 p t . The delta function in Eq. ( |4.2| ) tells that at a 
fixed number, u, of scatterings we have ((p 2 )) u = (((J2i=i k«) 2 )) = ((X^=i ^i))- The last 
equality is due to the azimuthal symmetry of the differential parton-parton cross-sections 
^jp. The resulting expression is symmetric under exchanges of the transverse momenta 
ki, so that ({p 2 ))u = u {{k 2 ))- Then it is immediate to see that 

(p 2 (x, b)) A = ±f d 2 p t dx'pAxx')T A (x>, b) = (p 2 (x, &)>i%^4r , (4.3) 
V A J d 2 p t Pa[x,o) 



where 



Jd 2 Pt dx'p 2d -^T A (x',b) 
^ X > 0))1 Jdx'a(xx')T A (x',b) 
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is the average transverse momentum squared in a single parton-parton collision. The pt- 
broadening of the wounded partons in a hA collision is then given by the p r broadening in 
a single collision multiplied by the average number of rescatterings suffered by a wounded 
parton. A similar result for the p t -broadening of a fast parton traversing a nuclear medium 
was derived in |27| . The p t - broadening was also studied in different contexts in [[], |T(| . 
Two interesting limits can be considered: 



(p?(ar,&)), 



(p 2 (x,b)) 1 aspo^oo 
{p 2 t (x,b)) 1 {n A (x,b)) aspo^O. 



(4.4) 



Since the minijet yield is dominated by transverse momenta of the order of the cutoff, 
these two limits say roughly that the minijets at high p t (i.e., high po in Eq. (|4.4|) ) suffer 
mainly one scattering. On the contrary, at low p t (i.e., low po in Eq. (|4.4| )) they undergo 
a random walk in the transverse momentum plane and the broadening is proportional to 
the average number of steps in the random walk, i.e., the average number of semi-hard 
scatterings suffered along the wounded parton trajectory. This picture will be studied in 



more detail in Sec. 5.1 



An explicit formula for the transverse spectrum can be obtained by studying its Fourier 
transform, since all the convolutions in Eq. ( |4.2|) turn into products and the sum over v 
may be explicitly performed. To this purpose, we introduce the Fourier transform of the 
parton-parton scattering cross-section 



o(y\ xx') 



d'k^^xx') 



Note that a(0;xx') = a(xx') and that due to the azimuthal symmetry of da/d 2 k, its 
Fourier transform depends only on the modulus, v , of v. Then, the transverse spectrum 



( |4.2| ) may be written as 

dWh 
d?p t 



{x,b,p) = r h (x,b-p) 



d 2 v 
(2nf 



-i\>t-vWh{v\x,i 



(4.5) 



where 



~ °° 1 T f 

W h (v;x,b) = / 

u=l V ' <-J 



dx'T a(x' ', b)a(v; xx' 



J dx'V A (x' ,b){5(v;xx')-&(0;xx')} 



- J dxT A (x',b)a(0;xx') 
-/ dx'r A (x',b)a(0;xx') 



(4.6) 



An immediate consequence is that the transverse spectrum has a finite limit as p t — >0, 
even when a cutoff on the momentum exchange is used: 



dW h 
d 2 p t 



Pt=0 



(x,b,P) = T h (x,b-P) 



d 2 v 
W) 



; W h (v;x,b) 
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4.2 Expansion in the number of scatterings 



We can obtain an expansion of Wh in the number of the rescatterings suffered by the 
incoming parton by expanding Eq. (|4.6|) in powers of a: 

oo 

W h (v;x,b) =J2Wh\v,x,b) 

— — \( I dxT^x', b) [a(v ; xx') — <x(0; xx')~\ j — ( — / dx'T a{x' ,b)a(Q]xx) 



(4.7) 



Coming back to the pt space, the expansion of the transverse spectrum in number of 
scatterings reads: 



The series Eq. (|4.7| ) can be obtained also by expanding W(v) around v = 0. Since the 
variable v is Fourier-conjugated to p t , the expansion of the transverse spectrum, Eq. (O), 



will be valid at high p t and we expect a breakdown of any truncation at sufficiently low 
momentum. Note that we can obtain this high-pt expansion of the spectrum directly in 
Pt space by expanding the exponential in ( f4.2|) and collecting the terms of the same order 
in a. As an example, the first three terms, Eqs. ( |A.1| ), ( |A.2| ) and ( |A.6| ), can be found in 



the appendix. The first two, suitably symmetrized, were used in 0, |3|, || || to explain the 
Cronin effect up to y/s = 39 GeV/A. The study of this series is the subject of Sec. 



numerical results up to n = 3 scatterings will be discussed in Sec. [5J] and compared to 
the whole spectrum. In the appendix we will discuss the symmetrization of the terms of 
the series. 



4.3 Cancellation of the divergences 

All terms of the expansion ( }4.8|) are divergent in the infrared region so that we need to cure 
them with the regulator po. Nevertheless, the infrared divergences are already regularized 
to a large extent by the subtraction terms originated by the expansion of exp[— (tia(x, b))} 
appearing in Eq. (|4.2|) , namely by the constraint of probability conservation. This cancel- 



lation mechanism was observed also in Ref . for the two-scattering term and in Ref . [|13 
in a different context. 

It is instructive to look in detail how the subtraction works for the lower order terms 
of the expansion. We start by considering the case of a single rescattering (u = 2). To 
simplify the notation we write the elementary differential cross-section da/d 2 k as cr(k), 
and notice that it depends only on the modulus, k, of the momentum. By expressing the 
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semi-hard cross-section as a = J d 2 ka(k) the term of order a 2 may be written as 

,(2) 



x, b, (3) = F h (x, b -P)f t a(x[, b)F A (x' 2 , b)dx' 1 dx' 2 d 2 k 1 d 2 k 2 (4.9) 

^(kx + k 2 - - 5^(1^ - Pt ) - <5< 2 >(k 2 - p t ) 



d 2 p 

<j(ki)o-(k 2 



x 

where the first term in the square brackets represent two successive scatterings with no 
absorption. The two negative terms are the corrections induced by the expansion of 
the absorption factor exp[— (n^(x, b))] of the single-scattering term, v = 1 in ( |4.2| ), and 
correspond to a single-scattering along with the effects of absorption in the initial or final 
state. The expression we obtained is symmetric in the integration variables ki and k 2 . 
The cutoff dependence is originated by the singular behavior of the integrand for ki pa 
or for k 2 ~ 0, since the 5-functions in the square brackets prevent the possibility of ki 
and k 2 being both zero at the same time. Because of the symmetry under the exchange 
ki <-> k 2 , to study the cutoff dependence of Eq. ( |4.9| ) it is enough to discuss the integration 
around k x = 0. In the region k x pa the term S^(ki — p t ) does not contribute, as long as 
p t is finite. The integration in k 2 is done with the help of the ^-functions and one obtains 



/ 



d 2 ha(k 



a(p t - ki) - o-(pt) 



On the other hand, for ki pa 0, one may use the expansion 

Pt ■ ki 



cr(p t - ki) ~ cr(pt) - a'(pt) 



Pt 



where p t ■ ki represents the scalar product of the two vectors, and cr'(p t ) = ^f^|°"(Pt) 
depends only on the modulus of p t . One is left with 



Pt 



p t ■ k^k^h = 



where the vanishing result is due to the azimuthal symmetry of cr(ki). The dominant 
contribution to the integral comes therefore from the next term in the expansion of 
a(ki — p t ), which goes as k\. Hence the resulting singularity is only logarithmic in 
Po, since cr(k) ~ k~ A as k — > 0. The subtraction terms, originated by the absorption fac- 
tor exp[— (ua(x, b))} in Eq. ( |4.9|) , have cancelled the singularity of the rescattering term 
almost completely. This feature is common to all the terms of the expansion ( [4.9|) . 

Hereafter we consider in detail the term with two rescatterings: 



d 2 p t 



(3) 



x, b, (3) = b — 13) J T yi (a; / 1 , 6)r^(x' 2 , 6)1^(23, b)dx' 1 dx' 2 dx' 3 d 2 kid 2 k2d 2 k 3 

x g (k 1 Mk 2 ) ff (k 3 ) r , 2)(ki + k2 + k ) 
o L 

- 5^(k! + k 2 - pt) - 5^ 2 \k 2 + k 3 - p t ) - 5 (2) (k 3 + ^ - p t ) 

+ 5^(k: - pt) + 5^ 2 \k 2 - P t) + 5^(k 3 - pt)} . (4.10) 
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The different 5-functions in Eq. Q4.10Q correspond to all the terms of order cr 3 in Eq. (|4.2|) 
and represent the triple scattering term together with all subtraction terms, induced by 
the expansion of the absorption factor exp[— (ua(x, b))] of the double- and of the single- 
scattering terms. The expression has been symmetrized with respect to k 1; k 2 and k 3 
and is singular for k x = 0, k 2 = and k 3 = 0. The 5-functions in ( 4.10|) prevent the 
tree momenta to be close to zero at the same time, then we start by discussing the most 
singular configuration corresponding to two integration variables both close to zero. Given 
the symmetry of the integrand it is enough to study the integration region with ki « 0, 
k 2 ~ 0. In this region the terms 5^(k\ + k 2 — pt), 5^(ki — p t ) and <5^(k 2 — p t ) do not 
contribute. The integrals on the transverse momenta are therefore written as 

• (4-11) 



d 2 k 1 d 2 k 2 a(k 1 )a(k 2 ) cr(p t - k x - k 2 ) - a(p t - ki) - cr(p t - k 2 ) + cr(p t ) 
In the region k x w 0, k 2 « one may use the expansion 



o- Pt - k -O-(Pt) - cr p t + - 

Pt 2 



ff "(P«) (p, ' k) 



Pt 



Pt 



(4.12) 



where p t X k represents the vector product of p t and k and cr"(p t ) = d \p t \s &('Pt) depends 
only on the modulus of pt- All terms proportional to cr(pt) cancel and all the terms linear 
in k integrate to zero thanks to the azimuthal symmetry of er(k). Then one is left with 



d 2 k\ d 2 k 2 cr ( kx ) o ( k 2 



2p 2 t 



2p 3 t 



(p t -(k 1 + k 2 )) 2 -(p t -k 1 ) 2 -( Pt -k 2 ) 2 " 
(pt x (ki + k 2 )) 2 - (p t x k x ) 2 - (p t x k 2 )" 



which simplifies to 

d 2 k 1 d 2 k 2 a(k 1 )cr(k 2 ) 



a " iPt) {pt ■ kj (p t ■ k 2 ) - ^(p t X k x ) ( Pt Xk 2 )|=0. 



Pt 



Pt 



The result is again zero because of the azimuthal symmetry of cr(k). Hence, all terms 



of the expansion ( 4.12Q up to the second order in k do not contribute. All other terms 
linear in ki or in k 2 , which are obtained from the first terms in the square brackets in 
Eq. (|4.11 ), do not contribute for the same reason, so the first term different from zero is 
at least of order kfk 2 , and originates a square-logarithm singularity as a function of the 
regulator p . 

One may repeat the argument for the regions where only one of the integration vari- 
ables is close to zero. We consider in detail the case kx w and k 2 , k 3 both finite. In this 
region the term 5^(p t — k) does not contribute to Eq. ( fl.lOp . The transverse momentum 
integrals are therefore 

J d 2 /cirf 2 A; 2 cr(k 1 )o-(k 2 )|cr(pt - kx - k 2 ) - a(p t - k x ) - a(p t - k 2 ) + cr(p t ) J 

+ J d 2 k 1 d 2 k 3 cr(k 1 )cr(k 3 )^-cr(p t - k x ) +a(p t )} . 
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To study the singularity it is sufficient to keep the first two terms in the expansion of 
cr(k) in |4.12| . One obtains 



Pt 



d 2 k\ d 2 k-2 a ( ki ) a ( k 2 ) 

v / i 1 ^ 'i , n (Pt - k 2) • ki p t ■ kx 

X S °"(Pi - K 2 ) - <7 (Pi - k 2 ) : cr(pt) + (T (p t ) CT{Pt ~ k 2 ) + cr( 

[ Pt - k 2 Pi 

+ | d^ ld 2^ (kiV(k3) |_ a(pt)+a , (pt) Pi^_ a(pt) | ) 

which simplifies to 

/ d^Aa^a^A-a'ivt - k 2 ) (Pt ~ ka) ' kl + 2a'{ Vt )^-\ = ■ 
J { Pt - k 2 Pt J 

As in the previous case one obtains a vanishing result thanks to the azimuthal symmetry 
of cr(k). All integrations in the singular points induce therefore at most a square-logarithm 
singularity, as a function of the cutoff, in the term with v = 3 in Eq. ( f4.8|) . 

The argument holds for the whole spectrum, as one may see by looking at its Fourier 
transform, Eq. ( f4.6|) . Indeed, to study the dependence of the inclusive spectrum on the 
regulator p at a given p t different from zero one needs to consider the first term in the 
square brackets only. The cutoff enters in the difference 



f da 


e ik-v _ x 


1 d?k 





2 2 tt f°° o da 

d k = —v — / k —T—dk + finite terms 
2 J po d*k 



so that, also in this case, the divergence for p — ► is only logarithmic. 



5 Numerical results and discussion 

In this section we discuss in detail, both qualitatively and quantitatively, the modifications 
induced by the rescatterings on the minijet inclusive transverse spectrum. We consider 
a proton-lead collision with center of mass energy ^Js = 6 TeV/A and impact parameter 
(3 = 0. In the numerical computations we used the leading order perturbative parton- 
parton cross-section with a mass regulator m = p Q : 

^;(-') = *pT§^«(-' s -V + m 2 )) 

where k is the fc-factor that simulates next-to-leading order corrections (we chose k — 2). 
The single-parton nuclear distribution function has been taken to be factorized in x and 
b: 

T A (x,b)=r A (b)G(x,Q) 
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where ta is the nuclear thickness function and G is the proton distribution function. We 
evaluated the strong coupling constant and the nuclear distribution functions at a fixed 
scale Q = m. In the computations we used a hard-sphere geometry 

r A (b) = A^VW^V6{R 2 - b 2 ) 

where R = 1.12A 1 / 3 is the nuclear radius. For G we used the GRV98LO parameterization 
PS| . At low p t the spectrum is obtained by computing numerically the Fourier transform 



in Eq. ( |4.5| ), but at high p t the result begins to oscillate too much, and in that region 
the spectrum was computed by using the expansion in the number of scattering up to 
the three-scattering term (the formulae actually used, Eqs. ( |A.1|) , ( |A.4|) and ( |A.8j ), are 



discussed in the appendix). We checked that the spectrum obtained by Fourier transfor- 
mation matched smoothly the expansion. 



5.1 Effects of rescatterings 

In this section we discuss the projectile and the target transverse spectrum averaged over 
a given rapidity interval: 

h (/3,rj min ,rj max ) = f dxd 2 b -^(x, 6, /3) , (5.1) 



« Pt Vmax ~ Vmin J n£[rj m i„,7)m.ax] Pt 

where we approximated the pseudo-rapidity by r\ = log(x-y/i/po)- The target spectrum, 
dWA/d 2 pt, is obtained by interchanging h and A in Eq. (|5 . 1| ) . Note that now we are taking 
into account all possible rescatterings of the target, as well. 

In fig.|I] we compare the full transverse spectrum (solid line) with its expansion in the 
number of scatterings up to three scatterings (dotted and dashed lines). We show both 
the projectile and target minijet spectrum in a pseudo-rapidity region rj 6 [3,4] for the 
projectile and rj 6 [—4, —3] for the target. Note that the rapidity is defined with reference 
to the projectile hadron direction of motion. The choice of a forward region (backward 
for the target) is done to enhance the effect of the rescatterings and to better discuss 
the deformation induced in the spectrum. Indeed, in those regions the average fractional 
momentum of an incoming parton is large, so that the number of available target partons 
is large and the probability of rescattering becomes large. 

First, we look at the projectile spectrum. At high p t the spectrum is enhanced with 
respect to the single scattering approximation because of the p t broadening induced by the 
rescatterings. As pt is further increased it approaches the single-scattering spectrum, as 
expected on general grounds when the pt distribution of the elementary scattering follows 
a power law. This can be understood qualitatively by looking at the path in p t -space 
followed by the incoming parton. Given a final large pt, due to the leading divergences in 
Eq. (|4.9| ), the leading processes to get that p t with two semi-hard scatterings are a first 
scattering with momentum transfer q\ « p Q followed by a second one with q 2 ~ Pt and 
vice- versa. For an analogous reason, the leading configuration to reach the final pt with 
three scatterings is q\ ~ p t plus g 2 ~ ?3 ~ Po and permutations. This sequence of three 
scatterings is less probable than the process with two scatterings as p t increases because 
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Pt ( Gev ) Pt ( Gev ) 

Figure 1: Left: Target pt-spectrum for rj E [—4,-3]. Right: Projectile p t -spectrum for rj 6 [3,4]. The 
full transverse spectrum (solid line) is compared with the one-, two- and three-scattering approximations 
(viz., dotted, short-dashed and long-dashed lines). 



the fraction of phase-space volume that this process occupies decreases much faster with 
p t than in the two-scattering case. For an analogous reason also the relative importance of 
the two-scattering term with respect to the single-scattering term decreases as pt increases. 
In conclusion as p t increases the average number of scatterings per parton decreases, and 
eventually the spectrum is well described by the single-scattering approximation. 

At intermediate p t the average number of scatterings per parton increases and the 
shape of the spectrum is more and more distorted with respect to the single-scattering 
case. In fact, the fraction of phase-space available to the leading configuration of a 
multiple scattering process (gi ss p t , qi ~ • • • ~ q n ~ Po and permutations) increases as p t 
decreases. However, this is not the only mechanism at work. Indeed, in our computation 
each wounded parton is counted as one minijet in the final state, independently of the 
number of rescatterings. On the other hand, in the single-scattering approximation one 
identifies the number of minijets in the final state with the number of parton-parton 
collision. This leads to an overestimate of the jet multiplicity and to a divergence of the 
spectrum at p t = as p goes to zero. Therefore at low p t the minijet yield is more and 
more suppressed with respect to the single scattering approximation. 

At very low transverse momentum p t < p a parton undergoes a large number of 
rescatterings, all with qi « p . Hence, the parton is doing a random-walk in the trans- 
verse plane and the spectrum becomes flat as p t — »0 because the phase space becomes 
isotropically populated. This shows that at very low pt multiple semi-hard scatterings are 



consistent with the random- walk model of Ref. flgl , while at moderate and high-j^ the 
physical picture is rather different. 

By comparing the results for the projectile and target transverse spectrum one sees 
that a projectile parton is traversing a very dense target and the effects of the rescatterings 
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p-Pb Vs = 6 TeV |J = fm m=l GeV K fact = 2 
: proj. : targ. : p + t : p + t (1 sc.) 




Pt ( GeV ) 

Figure 2: Projectile plus target pt-spectrum (solid line) at different rapidities compared to the result 
of the one-scattering approximation (dot-dashed line). Also shown are the contributions of the projectile 
minijets (dotted line) and and of the target minijets (dashed line). 



are large. On the contrary, a target parton sees a rather dilute system, and its minijet 
spectrum does not differ too much from the single-scattering result, except at very low 
p t . Moreover the changes induced by the rescatterings on integrated quantities, like 
those entering in the expression of the hadron-nucleus cross-section, are minimal. This is 
consistent with our approximation of not including rescatterings for the target partons to 
obtain analytic formulae for the hadron-nucleus cross-section. One can also see that the 
three-scattering approximation describes well the projectile spectrum for p t > 15 GeV, 
while it breaks down completely at p t < 7 GeV, where it becomes negative. For the target 
spectrum the three-scattering approximation is not accurate for p t < 4 GeV. 



5.2 Minijet inclusive transverse spectrum 

In this section we study the minijet transverse spectrum resulting from the sum of the 
transverse spectra of the projectile and target wounded partons: 

^ (Almm.u) = [ dxd 2 b (^-^(x,b,P) + ^-^-(x,b,P)) . (5.2) 



U Pt Vmax — Vmin Jri£[r] m i n ,r) max ] V « Pt « Pt 

We analyze the spectrum in three rapidity regions, namely rj e [—4,-3], 77 G [—1,1] 
and 77 G [3,4] (respectively "backward", "central" and "forward" with reference to the 
projectile direction of motion). While the target partons basically do not suffer any 
rescattering in all three regions, the projectile partons undergoes many rescatterings in 
the forward region, some in the central region and basically no one backwards. 

In Fig. we show the spectrum ( |5.2j ) (solid line) and the contributions of the projectile 
and of the target (dotted and dashed lines, respectively). For comparison also the total 
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Figure 3: Regulator dependence of the projectile plus target pt-spectrum at different rapidities for 
m= 1,2,3 GeV (viz., solid, dashed and dotted line). 

spectrum obtained in the one-scattering approximation is plotted (dot-dashed line). The 
spectra are computed with a regulating mass m = 1 GeV. 

In the backward region both the projectile and the target suffer mainly one scattering 
over all the p t -range except at p t ~ 0, and the spectrum is dominated almost everywhere 
by target minijets. 

In central and forward regions the target jets still suffer basically one scattering over 
all the p t range. On the contrary, the projectile crosses a denser and denser target and 
undergoes an average number of rescatterings that increases with pseudo-rapidity. This 
means that at low p t the projectile spectrum is very reduced with respect to the one- 
scattering approximation, and the minijet yield may become negligible with respect to 
the minijet yield from target. The overall effect is that at low p t the spectrum is dominated 
by minijet production from the target while at intermediate and high p t it is dominated 
by minijet production from the projectile. 

At very forward rapidities this effect becomes quite dramatic and the spectrum acquires 
a structured shape: it follows the inverse power behaviour of the single-scattering term 
at high p t , it is concave at intermediate p t because of the suppression of the projectile 
minijets and becomes convex again at low pt, where the target begins to dominate. 

In Fig.|3] we study the dependence of the spectrum on the choice of the cutoff, and 
plotted the result for m — 1, 2, 3 GeV. The deformation of the spectrum decreases as the 
regulator increases (indeed, the average number of rescattering decreases) and for m > 3 
GeV it begins to become negligible. 

The effects of the rescatterings are better displayed by studying the ratio of the full 
transverse spectrum and the single-scattering approximation: 




dW hA /d 2 p t dW hA /d 2 p t 



(5.3) 
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Figure 4: Ratio of the full projectile plus parton pt-spectrum to the one-scattering approximation at 
different rapidities and for m = 1, 2, 3 GeV (viz., solid, dashed and dotted line). 



where Ap ~ J d 2 bTh(b — /9)r^(6) is the number of target nucleons interacting with the 
projectile at a given impact parameter. 

In Fig.|] we plotted the ratio Rp(p t ), which measures the Cronin effect for minijet pro- 
duction, computed with three different regulators m = 1,2,3 GeV. At m = 3 GeV the 
effect of the rescatterings is rather small in all the three rapidity intervals, except at very 
low p t , and doesn't affect the integrated quantities like the average number of minijets. 
As the regulating mass is decreased the rescatterings begin to show up, and lead to a big 
effect in the forward region. 

The ratio Rp(p t ) is carachterized by three quantities: the momentum p x where the 
Rp crosses 1, the momentum pm where it reaches the maximum and the height Rm of 
the maximum. 

The sensitivity of p x on the cutoff decreases as the pseudo-rapidity increases. Loosely 
speaking, when the average number of scatterings is high, as it is the case at pt — Px, 
the jets loose memory of po> which gives the order of magnitude of the typical momentum 
exchanged in each collision, pu shows a slightly larger sensitivity on the regulator, since 
it lies in a region where the average number of scatterings is smaller. This behaviour is 
very different from the conclusions drawn by considering only the expansion up to two 
scatterings, where both p x and pm are proportional to po 0. In fact, at low center of mass 
energies the two-scattering is a good approximation in all rapidity ranges, except may be 
very forward. However, it breaks down in any case at transverse momenta comparable 
to the regulator p . Therefore, while most of the spectrum is well described by the two- 
scattering approximation, the behaviour of p x and pu is not. 

On the other hand, the height of the peak is much more sensitive to the cutoff, since its 
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leading term is roughly proportional to some power of the logarithm of the regulator: 



Since pu is not very large, the average number of rescatterings at that value of the 
transverse momentum, (n resc (pM)) , is much greater than one and the sensitivity of Rm 
on po is high. At high p t the average number of rescatterings tends to zero, so the 
sensitivity of the Rp on po decreases and disappears at very large transverse momenta. 
Note that the peak is located in a p^-region, where soft interactions (which have been 
disregarded in our approach) are expected to be negligible, therefore in that region our 
perturbative computations should describe almost completely the spectrum. Following 
Ref. || we might interpret po as the momentum scale at which the interaction deviates 
from the perturbative computations. With this interpretation p would acquire a physical 
meaning: though physics doesn't know about the artificial subdivision in hard and soft 
interactions, it is a well defined question to ask up to what scale are the perturbative 
computations good. If the collision dynamics would be determined by parton multiple 
elastic scatterings alone, then the measure of the height of the peak would be a way of 
measuring p . 

On the other hand, the sensitivity of Rp on p is rather signaling a lack in our description of 
the dynamics underlying the hadron-nucleus collision. We expect that such a sensitivity 
will be considerably reduced when including in the dynamics also the gluon radiation 
emitted by the multiply scattering partons. Some of the effects of the radiation on the 
transverse spectrum might be however described by the parameter p in the model where 
the radiation is neglected. Since the inclusion of gluon radiation in the dynamics would 
introduce new physical scales, like the radiation formation time, related to the energy of 
the collision and the nuclear size, we would expect in any case that the value of po will 
depend on y/s and A. 

6 Conclusions 

The purpose of the present article is to draw the attention to some of the advantages 
of studying hadron-nucleus semi-hard interactions at the LHC As in the case of lower 
energies, hA interactions represent an important intermediate step to relate hh and AA 
reactions, being much simpler to understand as compared with the latter. Moreover, 
even at higher energies, like those obtainable at RHIC and LHC, in hA collisions we 
don't expect the formation of a dense and hot system, like the quark-gluon plasma, so 
that one can study directly the nuclear modification of the dynamics without the need 
of disentangling the effects of the structure of the target and those due to the formation 
and evolution of the dense system. Hadron-nucleus interactions represent therefore the 
baseline for the detection and the study of the new phenomena peculiar to AA collisions. 

We faced the problem of unitarity corrections to the semi-hard cross-section by in- 
cluding explicitly semi-hard parton rescatterings in the collision dynamics, exploiting the 
self-shadowing property of the semi-hard interactions. In the interaction mechanism we 
took into account just elastic parton-parton collisions, while we neglected the production 



dW hA 
d 2 p t 





-J Pt=PM 
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processes at the partonic level (e.g., all 2 — > 3 etc. elementary partonic processes), whose 
inclusion represents a non-trivial step in our approach and deserves further study. 
Contrary to the case of AA collisions, we have been able to obtain closed analytic expres- 
sions for the semi-hard hA cross-section, Eq. (|3.7|) . To that purpose a crucial assumption 
has been to consider the hadron as a dilute system, so that rescatterings of nuclear par- 
tons can be neglected, while rescatterings of the projectile are fully taken into account. 
In our expressions we have disregarded correlations in the nuclear multi-parton distri- 
butions, whose effect may be nevertheless studied in a straightforward way within the 
present functional approach. 

We have then focused on the inclusive minijet transverse spectrum at fixed impact pa- 
rameter, Eq. ( |4.5|) , which is influenced in a more direct way by the rescatterings. The 
modifications of the transverse spectrum induced by the semi-hard rescatterings of the 
projectile partons is emphasized in the ratio Rp(p t ), Eq.|5T3|, defined as our p t spectrum 
divided by the impulse approximation. In particular, we have evaluated it at (5 — for 
different values of the regulator p . The results are described by the values of p x (defined 
by Rp{p x ) = 1), pm (which is the value of p t that maximizes the ratio) and Rm (which 
is the maximum of Rp). We obtain the both p x and pu depend weakly on p , while Rm 
has, on the contrary, a strong dependence on po also when the regulator is rather small. 
Therefore, the results for the spectrum allow also to identify the limits of the picture 
of the dynamics considered in this paper. Analogously to the average transverse energy 
and the number of minijets in AA collisions [F2TJ, some of the features of Rp, like p x and 
Pm, show a tendency towards a limiting value at small p . All these quantities depend 
therefore only marginally on details of the dynamics which have not been taken into ac- 
count in the present approach. On the contrary, the limits of the simplified picture of the 
interaction show up in R M . Because of the strong dependence of Rm on p , in order to 
describe the spectrum one needs in fact to fix experimentally the value of po by measuring 
Rm- This feature might be not so unpleasant, because if one limits the analysis to the 
inclusive transverse spectrum of minijets in hA collisions, all the effects which are not 
taken into account in the interaction (like the gluon radiation in the elementary collision 
process) are summarized by the value of a single phenomenological parameter. However 
this feature will not hold any further if one had to evaluate more differential properties 
of the produced state, which can be properly discussed only after introducing explicitly 
further details in the description of the elementary interaction process. 
The experimental measure of the Cronin effect in minijet production in hA collisions 
would be therefore of major importance: it would allow to establish the correctness of 
the whole approach here described and it would represent the basis for a deeper insight 
in the semi hard interaction dynamics both for hA and AA collisions. 
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A Expansion in the number of collisions 



For the numerical computation of the high-p f expansion of the minijet spectrum in the 
number of scatterings suffered by a projectile parton it is convenient to implement the 
subtraction of the IR divergences directly in the integrand. In this way the Monte-Carlo 
integrations, which we use because of the high dimensionality of the phase space (in 
particular for three or more scatterings), work at their best. In fact, Eqs. (|4.9| ) and (|4.10|) 
are not suited for numerical implementation due to the delta functions. The basic property 
that allowed the cancellation of the divergence in the integrand was the symmetry under 
exchanges of the integration variables. Unfortunately after using the delta-functions to 
perform on of the integrals, one obtains in general non-symmetric expressions. 

The goal of this appendix is to study how to symmetrize each term of the expansion of 
the transverse spectrum. We will discuss them in detail up to the three-scattering term, 
but the techniques discussed can be applied also to the generic term in the expansion. 
For simplicity, we will use the following notation, already introduced in the main text: 

-(k) = ■ 
A.l One-scattering term 

The one-scattering term doesn't include any subtraction term, so that we don't need to 
symmetrize it. It is simply given by 



dWh {1 \ 
d?p t 



x, b, (3) = T h (x, b-(3)[ dx'T A (x', 6)«7(pt) , (A.l) 



and corresponds to the result one obtains by considering just disconnected parton colli- 
sions and neglecting parton rescatterings. It corresponds also to modeling the hadron- 
nucleus collision as a superposition of hadron-nucleus collisions. 



A.2 Two-scattering term 



The two-scattering term is given by Eq. ( |4.9| ), and we need to perform one integration over 
ki or over k2 to dispose of the 5-functions. By calling simply q the remaining integration 
variable we obtain 



(x,b,/3) = T h (x,b-(3) I T A (x' 1 ,b)T A (x' 2 ,b)dx' 1 dx' 2 



x 



d 2 q 



0-(qV(Pt - q) - 2tr(q)<r(p t ) 



(A.2) 



As discussed in Section j4~3| , the negative term in the expression above subtracts the 
leading inverse power divergence in the integrand leaving only a logarithmic divergence. 
However, the cancellation happens only after performing the integral over q, which may 



25 



be a difficult result to achieve numerically (actually this is not a problem for the two- 
scattering term, due to the low dimensionality of the integral, but becomes a big issue 
from three scatterings on). 

There are two divergences to be subtracted: one in q ~ and the other in q ~ p t , 
but the subtraction term is divergent just in q ~ 0, and the cancellation of the inverse 
power singularities is obtained only after performing the integration over q. To allow 
the numerical integration to do a better and faster job, we want that the divergences in 
the convolution term and in the subtraction term be cancelled directly in the integrand. 
This is obtained by symmetrizing the integrand with respect to an interchange of the two 
singularities in the convolution term. Let's introduce therefore an operator that performs 
the interchange of the two singularities: 



T : q -> p t - q 



so that 



T J d 2 g/(q)= J d 2 g/( Pt -q) . 

Note that the change of variables operated by T has unit Jacobian and that T 2 = I. Then, 
we define the symmetrized two-scattering term as 



dW A 



(2) 



dW, 



(2) 



J (2) — A 
d 2 p t sym d?p- 

where we introduced the symmetrization operator 

1 



2 (I + T). (A.3) 



§(2) 

The result is: 

(x,b,/3) =Th(x,b — (3) / Ta(xi, b)YA{x' 2 , b)dx' t dx 

sym J 

x / d 2 q (i(q)a(pt - q) - er(q)a(p t ) - cr(p t - q)<r(p t ) . (A.4) 



dwf 



d 2 p t 



Note that the first term in ( |A.4|) describes two subsequent scatterings with total transverse 
momentum p t and is the naive pQCD result. The two negative terms are the absorption 
terms induced by probability conservation. The two IR divergences of the first term are 
canceled by these two subtraction terms: as q— >0 by the first one and as q— >p t by the 
second one. The remaining linear singularity gives a zero contribution because it is odd 
in a neighborhood of q = and q = pt so that only the logarithmic divergence remain. 
Note that now the two divergences are subtracted directly in the integrand, which was 
the goal of the symmetrization procedure. 

Eq. (|A.4j ) is the expression that we use in the numerical computations of the transverse 
spectrum at high p t . It could have been guessed directly from Eq. (|A.2|) , but the use of 
the symmetrization operator ( |A.3[ ) will facilitate the discussion of the more complicated 



three scattering term. 
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A. 3 Three-scattering term 



To prepare the ground for the treatment of the three-scattering term, we note that T 
generates the group of the permutations of the two singularities q ~ and q ~ p t ; this 
is called the symmetric group of order 2 and indicated as S 2 = (T) = {I, T}, where (T) 
means "generated by T". It's then easy to see that we can construct the symmetrizing 
operator ( |A.3|) by summing all the elements of 5*2 and by dividing by its cardinality. 

From Eq. ( 4.1U| ), after exploiting the (^-functions, the three-scattering term reads 



dW^ f 
„ h (x, b, (3) = T h (x, b - (3) I T A {x'ii b)T A (x' 2 , b)T A (x' 3 , b)dx' x dx' 2 dx'^ (A. 5) 

d% J " " 

x ^ J d 2 gd 2 r [a(q)a(r)a(p t - q - r) - 3a(q)a(p t - q)cr(p t ) + 3a(q)a(r)ff(p t )] . 

(A.6) 

Following the general analysis previously done at the end of the last paragraph, we observe 
that in ( |A.6| ) in absence of the cutoff we would have four divergences, i.e: 

q~0, r ~ 0, p t — q — r~0, p t — q ~ (A. 7) 

Then, to write the symmetrized three-scattering term, we need to consider the group 15*4 
of the permutations of these four divergences, which has 4! = 24 elements: 

<p(3) s(3)-p(2) 
' Bsym ~ ^ ' B 

where 

§(3) = iy T 

24 ^ 

TeS 4 

When applying this operator to the three-scattering term the resulting expression has 
49 terms and is too long to be discussed here. To have an idea of the result, we will 
consider only the subgroup S3 given by the permutations of the first three divergences 
in ( |A.71 ), which are the divergences that appear in the first term of (|A.6|) , i.e. the naive 



three-scattering term. After the symmetrization it will be immediate to check that all the 
"single" divergences cancel explicitly in the integrand, while "double" divergences cancel 
only after performing the integrations over the transverse momenta. We call "single" 
divergence a point (q, r) such that only one of the expressions in ( |A.7| ) is near zero, and 
"double" divergence a point such that two of these terms are nearly zero. For example 
{q ~ 0; r ^ 0,p t ,pt — q} and {q ~ 0; r ~ p t } are respectively a single and a double 
divergence. 

The first step is the definition of the operators that exchange the three singularities: 



Ti : < T 2 : <^ T. 






Note that they are idempotent: Tj = I. Next, we observe that the group S3 of the 
permutations of the three singularities is made of 3! = 6 objects, and that 

S 3 = (Ti,T 2 ,T 3 > = {To,T 1 ,T 2 ,T 3 ,T 4 ,T 5 } , 
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where T = I, T 4 = TiT 2 and T 5 = T1T3, so that the reduced symmetrizing operator is 



s(3) 



1 5 

3! ^ 

i=0 



Finally one can write the partially symmetrized three-scattering probability: 



dW 



(3) 



d 2 p t 



(x,b,P) 

sym 

T h (x, b - (3) 
1 



5 (3) - 

red 



dW 



(3) 



d?p t 



{x,b,P) 



Va(xi, b)TA(x' 2 , b)TA{x' 3 , fydx^dx^dx^d kid k 2 d k. 



3! 



x — / dx'd 2 qd z r 



tr(q)<r(r)<r(pt - q - r) 



~a(q)a(r)a(p t - q) + -a(q)a(p t - q)tr(p t ) 



-a{q)a(r)a(p t 



1 f 



<r(pt - q - r)a(q)cr(p t - q) + -<r(p t - r)<r(p t - q)cr(pt) 



0"(Pt - q - r)tr(r)tr(p t - r 

q 



-o-^rjo-^pt 



r)a{q + r) + -a{r)a{q_ 



r)a(p t - q 
- r)o-(pt) 



-(T(q)o-(pt - q - r)cr(q + r) + -cr(q)(r(q - r)cr(p t ^ 



(A.8) 



Analogously to what has been done for the two-scattering term, one can see by inspection 
that the four single divergences ( |A.7| ) explicitly cancel in the integrand, while double 
divergences cancel only after performing the integrations over q and r. By considering all 
four singularities, and by using the whole 64 group we would get explicit cancellation of 
both "single" and "double" divergences directly in the integrand. Nonetheless, the partial 
symmetrization is enough to get satisfactory numerical results. 

In conclusion, to compute numerically the expansion of the transverse minijet spectrum 
in the number of scatterings one has to fully exploit the symmetry properties of each term, 
in such a way that all the divergences get cancelled directly in the integrand. This is crucial 
to obtain a good numerical precision and to speed up the computation of the terms with 
three or more scatterings. In this appendix we developed a general technique to perform 
such a symmetrization. 
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